Reentrant superconductivity in conical-ferromagnet/superconductor nanostructures 
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We study a bilayer consisting of an ordinary superconductor and a magnet with a spiral magnetic structure 
of the Ho type. We use a self consistent solution of the Bogolioubov-de Gennes equations to evaluate the pair 
amplitude, the transition temperature, and the thermodynamic functions, namely, the free energy and entropy. 
We find that for a range of thicknesses of the magnetic layer the superconductivity is reentrant with temperature 
T: as one lowers T the system turns superconducting, and when T is further lowered it turns normal again. This 
behavior is reflected in the condensation free energy and the pair potential, which vanish both above the upper 
transition and below the lower one. The transition is strictly reentrant: the low and high temperature phases are 
the same. The entropy further reveals a range of temperatures where the superconducting state is less ordered 
than the normal one. 



More than thirty years ago, reentrant superconductivity as- 
sociated with magnetic ordering was first observed in the 
ternary rare-earth compounds ErRli4B4 and HoMogSgfl— 5]. 
On cooling, these materials first become superconducting at 
a critical temperature Upon further cooling, inhomoge- 
neous magnetic order sets in. This ordering coexists with 
superconductivity [6] over a very narrow T range. This on- 
set is nearly immediately [7] followed by that of long-range 
ferromagnetic order, which entails the destruction of super- 
conductivity, at a second critical temperature T c \. Thus, the 
reason for the disappearance of the superconductivity at T c \ 
is essentially the presence of the magnetism. That nonuni- 
form magnetic ordering can appear in the presence of super- 
conductivity is consistent with the prediction made by Ander- 
son and Suhl[8]. Reentrant superconductivity of a different 
kind is also found in ferromagnet/superconductor (F/S) lay- 
ered heterostructures[6]. On increasing the thickness, dp, of 
the ferromagnet layers in such structures, while keeping the 
thickness of the superconductor layers constant, the super- 
conductivity may disappear for a certain range of thickness 
(dpi < dp < dfi) and then return for larger dp (dp > dpi). 

The purpose of this Letter is to show that superconductivity 
in F/S nanostructures which is reentrant with temperature can 
occur under some circumstances, when the magnetic structure 
is non-uniform. That is, for certain types of ferromagnets, the 
Cooper pair amplitude in such structures can be nonvanishing 
in a range T c i < T < T C 2, with T c \ > 0. Specifically, we 
have found that this reentrance occurs in F/S bilayers where 
the magnetic order of the F layer is of the spiral type, as in 
Holmium[9]. The reentrance we find is very different from 
that in ErRh4B4 or HoMoeSs. There, the high T phase is para- 
magnetic and the low T phase is ferromagnetic. In our case, 
the magnetic order remains unchanged: it is the same above 
T c \, below T C 2, and in between. Reentrance occurs also[10] 
in some quasi one dimensional superconductors, but there the 
low T phase is insulating. In our case, we have strict reen- 
trance: the lowest T and highest T phases are the same, while 
in the entire range in between, superconductivity and mag- 
netism harmoniously coexist. This is unusual. Superconduct- 
ing reentrance is also found in granular films[ll]: it is not 



due to magnetism but it involves the turning on and off of the 
intergrain Josephson coupling. Here, we are able to evaluate 
the thermodynamic functions of the system as it undergoes 
the transitions, and from their behavior one can glimpse the 
reasons for the occurrence of the reentrance. The balance be- 
tween the internal energy of the system and its entropy can 
result in a situation where the entropy of the thermodynam- 
ically stable superconducting state is higher than that of the 
normal state. 

Extensive theoretical[6, 12-17] work indicates that the ori- 
gin of dp reentrant superconductivity in F/S nanostructures 
can be traced to the damped oscillatory nature of the Cooper 
pair wave functions in ferromagnets [18, 19]. Qualitatively, 
when a Cooper pair enters into an F region, it decays and 
the electron with magnetic moment parallel to the internal ex- 
change field h lowers its energy by an amount proportional to 
h, while the other electron with opposite spin raises its energy 
by the same amount. Then, the kinetic energy of each elec- 
tron changes and as a result[18] the Cooper pair entering into 
an F region acquires a spatially dependent phase in the F layer. 
This propagating character of the Copper pair leads to interfer- 
ence between the transmitted pairing wave function through 
the F/S interface and the reflected wave from the opposite sur- 
face of the ferromagnet. Experimentally, the reentrant behav- 
ior of superconductivity with dp has been observed and con- 
firmed in Nb/Cui_, v Ni A bilayers and Fe/V/Fe trilayers [20-22]. 
However, in the work we present here, reentrance occurs with 
temperature, rather than just with geometry. Thus, although 
it is already known[23] that the nonuniform Ho structure has 
strong effects in the S/F proximity phenomena, no T reen- 
trance results have been predicted or observed. 

In the rest of this paper, we will first review our methods 
as applied to Holmium/superconductor (Ho/S) structures and 
then discuss the microscopic behavior of the pair amplitudes 
as well as the thermodynamic quantities. The approach we 
use here is based on exact, self-consistent, diagonalization of 
the Bogoliubov-de Gennes (BdG)[24] equations for clean F/S 
structures. This approach not only has the virtue of being very 
general but is also able to describe short wavelength oscilla- 
tions, which is important for small structures. The self consis- 
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FIG. 1. (Color online) Schematic of the ferromagnet (Ho) - super- 
conductor (S) bilayer studied. The conical ferromagnet has a spiral 
magnetic structure described by an exchange field h, (see text). The 
system is infinite in the x — z plane and y is normal to the interfaces. 



tent methods we use to diagonalize the BdG equations have 
been extensively described in the literature (see e.g. Ref. 25 
and references therein) and details will not be given here, ex- 
cept where crucial. 

The geometry of the Ho/S system we consider is depicted 
schematically in Fig. 1. The y axis is normal to the lay- 
ers. The system is assumed to be infinite in the x-z plane 
and has a total length d in the y direction. The S layer in 
our assumed Ho/S system is a conventional s-wave supercon- 
ductor with thicknesses ds and a Ho layer of thickness dF- 
As in previous work, the magnetic structure is described via 
a local exchange field h which in this case is of the form: 
h = Iiq Jcos By + sin [sin ( ^ J x + cos (&\ z] J, where for Ho 
we have[9, 23] = An/9 and tp = tt/6. We will take a, the 
lattice constant, as our unit of length and assume throughout 
that the system is below the temperature (21 K) at which, 6 
switches from n/2 to An/9, i.e. Ho becomes ferromagnetic. 

The effective Hamiltonian, r H e ff< that we use to model our 
Ho/S system takes the form 

Keff = J d 3 r[J] ^(r) {-^ - £/) -Mr) 
-^^(rKh-^r)}, 



(1) 



where A(r) is the usual singlet pair potential; t//' a and \p a are 
the creation and annihilation operators with spin a respec- 
tively; Ef is the Fermi energy and er are the Pauli matrices. 
To recast the < H e ff into diagonal form, we apply a generalized 
Bogoliubov transformation, if/ a (r) = £ [M„Q.(r)y„ + v* [t (r)y]J, 

n 

where the quantum number n enumerates the quasiparticle 



(u„ a ) and quasihole (v mr ) spinors. The y„ and yj, are the 
Bogoliubov quasiparticle annihilation and creation operators 
respectively. By making use of the commutation relations, 
1~(eff,y„] = -e n y„ and \^H e ff,ji\ = e„y„, one obtains the 
3dG equations in matrix form. In the geometry chosen, the 
dependence of the wavefunctions on the x and z variables 
leads to an obvious phase factor that can be canceled out. This 
results in a set of quasi one dimensional problems of the form: 
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where H e = -(1 /2m*)(3 2 /dy 2 ) + e ± -E f , with e ± being the 
kinetic energy associated with the transverse direction. Thus 
the spatial dependence of the amplitudes is only on y. The 
exchange field h(y) in Ho is nonvanishing only in the F region 
andprecesses as given above (see also Fig. 1). The pair poten- 
tial must be determined self-consistently by solving the BdG 
equations together with the condition, 



AGO = ^J)' - u nl(y»*„M tanh(^), (3) 
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where T is the temperature, and g(y) is the usual BCS cou- 
pling constant associated with a contact potential that exists 
only in the S region. The prime on the sum implies that only 
states corresponding to positive energies below the "Debye" 
cutoff a>D are included. The self consistent diagonalization is 
achieved as in the previous work mentioned above, the only 
difference being that the matrices to be diagonalized are in 
this case unavoidably complex. 

From the self consistent results one can evaluate immedi- 
ately the pair amplitudes and, as explained below, the thermo- 
dynamic quantities. The transition temperatures can be most 
conveniently evaluated by a linearization method[26]. Near 
the transition temperature, the equation for A can be written 
as A, = Yjq JiqAq, where A, are the expansion coefficients with 
respect to the orthonormal basis, (pi(y) = y2/d sm(iny/d), and 
J iq is given as J iq = (/" + J\ ' )/2, where 



(4) 



(5) 



Here y = yo/4nD with yo being the dimensionless coupling 
constant in S; D is the total dimensionless thickness of the 
structure, D = k/sd, and kf$ is the Fermi wavevector in 
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FIG. 2. (Color online) Calculated transition temperature T c , nor- 
malized by T° (see text), vs. the dimensionless ferromagnet width, 
D F (= dfkfs ). Main plot: The upper points ((red) + , (green) x signs) 
are the usual critical temperature (T c2 ), leading to the superconduct- 
ing state as T is lowered. In the region 4 < D F < 5 (highlighted 
by the (green) x signs) a second transition back to the normal state 
appears at the (blue) star points forming the lower "dome". The in- 
set shows a broader range of magnet widths, revealing the overall 
periodicity of T c2 . 



S. We take kf S = I /a; e"^'° are unperturbed particle(hole) 
energies; and F inm = n V2d % pq (u° npt u° m ql - u^u ^) K irm , 
G inm = ^EwK/^-v^J^, where K inm = 

J dyg(y)(f>i(y)(p n (y)(/> m (y). The u° np and v° m are the expansion 
coefficients of the unperturbed (A = 0) particle (hole) ampli- 
tudes in terms of the basis set. 

This linearization method is easily used to evaluate the tran- 
sition temperature. As explained in Ref. 26, one simply has 
to find the largest eigenvalue, A, of the matrix J lq and see if it 
is greater or smaller than unity: in each case one is, respec- 
tively, in the superconducting or the normal state. The tran- 
sition temperatures are those at which the largest eigenvalue 
changes from greater to smaller than unity: one finds T c by 
evaluating A as a function of T. In the usual case A is smaller 
than unity when T is larger than T c . In a reentrant case with 
superconductivity in the range T c \ < T < T c o, we find T c \ by 
increasing T from zero until A > 1 and T C 2 by decreasing T 
from above T C 2 until A > 1 . 

In all results given here, the thickness of the S layer is fixed 
at ds = (3/2)^o, where £o is the usual BCS coherence length 
in S. We take £o = lOOA:^], and vary dp. The magnitude of h 
is Q.\5Ef. Results for the transition temperature, normalized 
to the bulk transition temperature T° of S, are shown in Fig. 2, 
plotted as a function of Df = dfk/s . In the inset, we see that 
the overall behavior of T c consists of the expected damped os- 
cillations with approximately the Df periodicity of the spiral 
magnetic structure (twelve, in our units). The main plot shows 
in more detail the structure near the first minimum. There we 
see also a lower small dome-shape plot ((blue) stars) with a 
maximum at Df as 4.5. The system is in the normal phase 



0.5 
0.45 
0.4 

«|? 0.35 
[T 0.3 



E 
o 



0.2 
0.15 
0.1 

0.05 




8 






5 o 

CO 


_y \ > 


. - 


« -4 

-a 








3 0.2 0.4 

T7T C ° 








- 




■ ■ 
■ ■ 









0.0002 



-0.0002 



- -0.0004 <j 



-0.0006 



-0.0008 



-0.001 



0.1 0.2 n 0.3 0.4 0.5 

T/T° 



FIG. 3. (Color online). Pair amplitude and thermodynamic func- 
tions. All quantities are plotted vs. T/T°. In the main plot, the 
(red) triangles and left vertical scale display the normalized (see 
text) singlet Cooper pair amplitude F(Y), one correlation length in- 
side S. This quantity vanishes at the upper transition temperature 
(T c2 ~ 0.47Tf) and again at the lower transition T cl « 0.07T®. The 
(blue) squares and right scale are the normalized (see text) conden- 
sation free energy, A/. The vanishing of A/ at the upper and lower 
transitions is clearly seen. The inset shows the normalized entropy 
difference A5 = -(dAf /d(T /T°)). 



inside the dome and, at constant Df, it is in the superconduct- 
ing phase between the two curves. In the Df range including 
the dome, the system, upon cooling, first becomes supercon- 
ducting at a higher temperature T c n, and with further cooling, 
returns to the normal phase at a lower temperature T c \. 

In Fig. 3 we display additional direct evidence confirming 
the existence of the reentrant behavior and showing its prop- 
erties. All results in the figure are for a system in the reentrant 
region, with Df = 4.3, and are plotted vs. T/T® . We con- 
sider first (main plot, (red) triangles, left vertical scale), the 
Cooper pair amplitude F(y) defined by A(y) = g(y)F(y) (see 
Eqn. (3)). The quantity shown is F(y = go), normalized to 
its bulk value in S, at a position one coherence length inside 
S. This amplitude vanishes below T c \ and above T C 2, with the 
values of T c \ and agreeing with those previously found: 
we can see from Fig. 2, T cl « 0.077; and T c2 ~ 0.47r° at 
Df =4.3. The continuity of the pair amplitude at T c \ and 
also indicates that the transitions are of second order. 

In the rest of Fig. 3 the thermodynamics of the transitions, 
which follows from the free energy, is shown. Using a stan- 
dard formalism [26, 27], we calculated Fs, the free energy of 
the whole system in the self consistent state, and Fn, the nor- 
mal state (A = 0) free energy. The normalized condensation 
free energy A/ = {Fs - Fn)/(2Eo) (Eq is the condensation 
energy of bulk S material at T — 0) is then plotted in the main 
part of Fig. 3 ((blue) squares, right scale). Both Fs and F^ are 
monotonic and have negative curvature with T as required by 
thermodynamics, but their difference is nonmonotonic. Al- 
though Af is small compared to its bulk value, we can still 
identify the two transition temperatures T c \ and T C 2 from this 
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plot. Their values are again in agreement, within numerical 
uncertainty, with those found from the pair amplitudes and 
from direct calculation. The system is in the superconduct- 
ing state when the T falls in the range T cl < T < T C 2- As 
r c i is approached from above or T C 2 from below, the solu- 
tion with A J disappears (as seen in the amplitude plot, 
(red) triangles), and the two free energies coincide: this is 
just what happens in ordinary BCS theory as the transition 
is approached from below. The minimum condensation free 
energy occurs at T m « 0.32T' which coincides with the lo- 
cation of the maximum pair amplitude. We also evaluated the 
entropy in the normal and superconducting states via textbook 
formulas. The normalized [26] entropy difference for the same 
case is shown in the inset of Fig. 3. It confirms that the system 
indeed undergoes second order phase transitions at both T c \ 
and T C 2. Unlike in a bulk superconductor, or in non-reentrant 
structures [26], there is now a range of T (T c \ < T < T m ) 
where the superconducting state is less ordered than the nor- 
mal one, and the entropy helps maintain the superconductivity. 

What is the physics behind this T reentrance? For F/S bi- 
layers with a uniform ferromagnet, the superconductivity dis- 
appears for a certain range of dp. This disappearance is due to 
the oscillating Cooper pair amplitude. Now, the spiral magne- 
tization in Ho introduces an oscillating magnetic order. Both 
the magnetic structure and the superconductivity are nonuni- 
form, consistent with the prediction in Ref. 8 that supercon- 
ductivity may coexist with nonuniform magnetic order. Ther- 
modynamically, we have here a subtle example of entropy- 
energy competition. In the range T m < T < T C 2, Af and AS 
behave qualitatively as they do for an ordinary [28] bulk super- 
conductor in the region < T < T c , (although they are much 
smaller). In either case AS vanishes at both ends of the range 
and has a minimum at a finite T in between. But in our case 
T m is nonzero. For T < T m , AS turns positive because of the 
oscillatory nature of the pair potential. The superconducting 
state becomes then the higher entropy phase: the roles of the 
N and S phases are thus reversed, the pair potential begins to 
decrease, and this leads inexorably to the lower transition at 
T c \, and to the reentrance into the same N phase. 

We have already seen above the clear differences between 
this entropy competition driven situation and other singlet su- 
perconducting T reentrance cases associated with field in- 
duced situations. Temperature reentrance involving long 
range magnetic order has been long known to occur in spin 
glasses [29], but the lowest T and high T phases (spin glass 
and paramagnetic respectively) are not the same. Somewhat 
similar but even more complicated situations occur in liquid 
crystals and may be a general property of[30] frustrated sys- 
tems. But a survey would take us too far afield. 

To our knowledge, this effect has not been searched for. 
The predicted range of T needed, down to about 0. 1 should 
pose no difficulty. The best course should be to fabricate sam- 
ples of varying dp, verify the T c oscillations (see Fig. 2 in- 
set) and then search for reentrance for dp near a minimum 
of the T c vs dp curve, where the phenomenon is predicted to 
occur. (This is possibly because such minima are associated 



with fragility of the superconducting state). It has proved ex- 
perimentally feasible[31] to study the T induced - n state 
transitions in S/F/S trilayers, which are related to a different 
effect[26] also involving nontrivial pairing correlations. Thus, 
that difficulties in sample making are not insurmountable. 

In conclusion, we predict that F/S bilayers with an inho- 
mogeneous conical magnetization will exhibit reentrant su- 
perconductivity with T, in addition to dp. Thus, superconduc- 
tivity exists for T c \ < T < T C 2 with nonzero T c \ under some 
conditions. We have shown clear evidence for this by self con- 
sistently determining the critical temperature-thickness phase 
diagram, and the T dependence of the pair amplitude. The 
thermodynamics were investigated via the free energy, reveal- 
ing a range of temperatures in which the normal state is lower 
in entropy than the superconducting state. 
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